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Summary 

Cell cycle events are driven by Cyclin dependent kinases 
(CDKs) and by their counter-acting phosphatases. Activation 
of the Cdkl: Cyclin B complex during mitotic entry is 
controlled by the Weel/Mytl inhibitory kinases and by 
Cdc25 activatory phosphatase, which are themselves 
regulated by Cdkl: Cyclin B within two positive circuits. 
Impairing these two feedbacks with chemical inhibitors 
induces a transient entry into M phase referred to as 
mitotic collapse. The pathology of mitotic collapse reveals 
that the positive circuits play a significant role in maintaining 
the M phase state. To better understand the function of these 
feedback loops during G2/M transition, we propose a simple 
model for mitotic entry in mammalian cells including spatial 
control over Greatwall kinase phosphorylation. After 
parameter calibration, the model is able to recapture the 
complex and non-intuitive molecular dynamics reported by 



Potapova et al. (Potapova et al., 2011). Moreover, it predicts 
the temporal patterns of other mitotic regulators which have 
not yet been experimentally tested and suggests a general 
design principle of cell cycle control: latching switches buffer 
the cellular stresses which accompany cell cycle processes to 
ensure that the transitions are smooth and robust. 

© 2013. Published by The Company of Biologists Ltd. This is an 
Open Access article distributed under the terms of the Creative 
Commons Attribution License (http://creativecommons.org/ 
licenses/by/3.0), which permits unrestricted use, distribution 
and reproduction in any medium provided that the original 
work is properly attributed. 
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Introduction 

The eukaryotic cell cycle is the process during which a mother 
cell duplicates all its components and divides them more or less 
evenly between two daughter cells, providing them with the 
necessary information and machinery to repeat the process. In 
order to keep the genetic content constant over subsequent 
generations, chromosomes are alternately replicated in S phase 
and then equally segregated between two daughter cells in M 
phase. Between these two cell cycle events are two gap phases 
(Gl and G2) before and after S, respectively. Both Gl and G2 
allow the cell to assess the readiness to proceed to the next cell 
cycle phase (Morgan, 2006). 

Cell cycle events are controlled by phosphorylation through 
Cyclin dependent kinases (CDKs; Cdkl, Cdk 2, Cdk 4 and Cdk6) 
and dephosphorylation catalysed by counter-acting phosphatases 
(e.g. PP2A-B55). CDKs require activatory subunits called 
Cyclins (Cyclin A, Cyclin B, Cyclin D, Cyclin E etc.), most of 
whose levels oscillate during the cycle. In humans, Cyclin B is 
the major activatory subunit required by Cdkl to promote mitosis 
and most of Cyclin B (CycB) is synthesised over G2 (Brizuela 
et al., 1989). However, during S and G2 phase the Cdkl subunit 
of the complex is phosphorylated (Draetta et al., 1988), which 
inhibits Cdkl: CycB activity and stops premature entry into 
mitosis (Krek and Nigg, 1991). Inhibitory phosphorylations are 
carried out by two kinases, Weel and Mytl. Weel preferentially 
phosphorylates and inhibits Cdkl at Tyrosine 15 (McGowan and 



Russell, 1993; Parker et al, 1995; Watanabe et al., 1995). 
However Weel is phosphorylated by Cdkl during mitosis, which 
inhibits its kinase activity (Parker et al., 1995; Watanabe et al., 
1995). Mytl preferentially phosphorylates and inhibits Cdkl on 
Threonine 14 and is also inhibited by Cdkl phosphorylation in 
mitosis (Booher et al., 1997). CdkhCycB also phosphorylates and 
activates an activatory phosphatase, called Cdc25, during mitotic 
progression. Cdc25 removes the inhibitory phosphorylations 
on Cdkl (Hoffmann et al., 1993) and its depletion can cause G2 
arrest in human cells (Galaktionov and Beach, 1991). Recent 
experiments suggest that M phase entry requires concomitant 
inhibition of a Cdkl counter-acting phosphatase (PP2A-B55), 
which is achieved by Greatwall-kinase dependent activation of 
phosphatase inhibitors (Gharbi-Ayachi et al., 2010; Mochida et al., 
2010). 

The feedbacks between Cdkl:CycB and the enzymes regulating 
Cdkl inhibitory phosphorylations (Weel/Mytl and Cdc25) control 
the activation of Cdkl:CycB. These feedbacks create switch-like 
activation and inactivation thresholds for Cdkl:CycB. The 
activation and inactivation thresholds of Cdkl:CycB for mitotic 
entry and exit are not the same, a phenomena called hysteresis 
(Novak and Tyson, 1993). The M phase state can be maintained 
with a lower Cdkl activity than the activation threshold, therefore 
M phase is a self-maintained state. Hence, mitotic entry is an 
irreversible process controlled by a bistable switch (Pomerening 
et al., 2003; Sha et al., 2003). Removal of the bistable switch in 
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human cells causes rapid cycles between a mitotic like state 
without proper exit and cytokinesis (Pomerening et al., 2008). 

In a recent paper, the functional significance of these 
positive feedbacks was directly tested in HeLa Cells. When 
both Weel/Mytl and Cdc25 were inhibited in early G2 phase, 
cells prematurely entered into mitosis. However, cells were 
unable to maintain their mitotic state after nuclear envelope 
breakdown (NEBD) and the mitotic phenotype was reverted. 
Compromising Cdkl's activatory positive feedback loops by 
inhibition of Weel/Mytl and Cdc25 makes mitotic entry 
reversible, and this phenomenon has been coined 'mitotic 
collapse' (Potapova et al., 2011). 

Mitotic collapse clearly indicates that cells experience some 
challenges during mitotic entry. Normally, the challenges are 
handled by the positive feedbacks through Cdc25 and Mytl/ 
Weel. However if the feedbacks are weakened, the cellular 
challenges could result in serious problems. The phenomenon of 
mitotic collapse raises three related questions. What could be the 
cellular challenges during mitotic entry? How do cells normally 
handle these challenges and why do cells fail to do so during 
mitotic collapse? 

To answer these questions, we present a dynamical model 
for the G2/M transition of mammalian cells. To handle the 
complexity brought about by NEBD, the model considers spatial 
control of Greatwall phosphorylation, which allows us to 
recapitulate the intriguing dynamics reported by Potapova et al. 
in quantitative details (Potapova et al., 2011). We propose that 
NEBD might cause a cellular stress during mitotic entry, which 
would usually be buffered by the bistability in Cdkl activation. 
When the bistability is compromised (such as by loss of positive 
feedbacks through Weel/Mytl and Cdc25), the cellular buffering 
capacity is impaired and mitotic entry becomes reversible and 
mitotic collapse takes place. This case-study reinforces the claim 
that bistability is an essential property for irreversible cell cycle 
transitions (Novak et al., 2007). 

Model construction 

The wiring diagram 

The molecular mechanism of our model used to analyse the 
kinetics of mitotic entry of HeLa cells in the experiments of 
Potapova et al. (Potapova et al., 201 1) is summarized in a wiring 
diagram (Fig. 1). In the experiments, nocodazole was used to 
disrupt the mitotic spindle which activates the mitotic checkpoint 
that blocks CycB degradation by Anaphase Promoting Complex/ 
Cyclosome (APC/C). To simplify the model, we assume that 
CycB is only synthesized in S/G2 phases and its level becomes 
constant once cells enter into mitosis. Since Cdkl is in excess of 
CycB (Arooz et al., 2000), the formation of Cdkl:CycB complex 
is limited by CycB availability, therefore we only monitor CycB 
levels which will correspond to the Cdkl:CycB complex. Mytl/ 
Weel kinases and Cdkl:CycB phosphorylate and inhibit each 
other in a double-negative feedback (Booher et al., 1997; Parker 
et al, 1995; Watanabe et al, 1995). Meanwhile, Cdc25 and 
Cdkl:CycB activate each other in a positive feedback loop 
(Hoffmann et al., 1993). 

CdkLCycB also phosphorylates and activates Greatwall kinase 
(Blake-Hodek et al., 2012; Yu et al., 2006), which is abbreviated 
as Gwl in this model. Greatwall-kinase phosphorylates and 
activates two phosphatase inhibitors, Endosulphine and Arppl9 
(Gharbi-Ayachi et al., 2010; Kim et al, 2012; Mochida et al., 
2010). These two inhibitors are referred here simply as ENS A. 
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Fig. 1. The wiring diagram of the mitotic control model. Weel/Mytl and 
Cdkl:CycB inhibit each other, while Cdc25 and Cdkl:CycB activate each 
other. Cdkl:CycB phosphorylates Greatwall kinase (Gwl). The phosphatase 
PP2A-B55 is assumed to dephosphorylate Weel/Mytl, Cdc25 and Gwl. The 
active Gwl-P activates ENSA and thus represses PP2A-B55. When ENSA is 
dephosphorylated, it no longer binds to PP2A-B55. A generic mitotic substrate 
(S), phosphorylated and dephosphorylated by Cdkl:CycB and PP2A-B55, 
respectively, is an indicator of the mitotic progression. 

The phosphorylated forms of these inhibitors (ENSA-P) bind and 
inhibit PP2A-B55, which dephosphorylates CdkLCycB target 
proteins, including Weel/Mytl and Cdc25 (Clarke et al., 1993; 
Mochida et al., 2009). 

Since Greatwall is a Cdkl:CycB substrate, it has been 
proposed that it is dephosphorylated by PP2A-B55 (Domingo- 
Sananes et al., 2011). Though this assumption has not been 
directly confirmed by experiments yet, it is supported by the 
following experimental evidence. If protein synthesis is repressed 
and Cdc25 is depleted in an interphase Xenopus extract, Cdkl 
activity is low and Gwl is unphosphorylated. Despite of the low 
Cdkl activity, Gwl becomes abruptly phosphorylated when 
PP2A-B55 is repressed by addition of phosphorylated 
Endosulfine or Arppl9 (see figure 2B of Gharbi-Ayachi et al. 
(Gharbi-Ayachi et al., 2010)). This experiment implies that 
PP2A-B55 dephosphorylates Gwl. Therefore Gwl and PP2A-B55 
inhibit each other which adds another positive feedback loop to 
the network. 

Spatial controls and the effect of nuclear envelope breakdown 
The majority of Weel is located in the nucleus (Heald et al., 
1993) but Mytl is located in the cytoplasm (Liu et al., 1997). In 
this regard, an inhibitory kinase is available to phosphorylate 
Cdkl in both the nucleus and cytoplasm. Hence, in the current 
model we use Weel/Mytl to represent the overall Cdkl 
inhibitory kinase activity. Similarly, Cdc25C is available in 
cytoplasm (Dalai et al., 1999) while Cdc25A is in the nucleus 
(Kallstrom et al., 2005), so we use Cdc25 to represent their 
cumulative phosphatase activity. 

CdkLCycB is activated and imported into the nucleus just 
before nuclear envelope breakdown (Gavet and Pines, 2010; 
Santos et al., 2012). However, significant CdkLCycB remains in 
the cytoplasm, which could be important to ensure 
synchronisation of cytoplasmic and nuclear mitotic events 
(Gavet and Pines, 2010). We assume that active CdkLCycB 
shuttles between the cytoplasm and the nucleus, being enriched in 
the nucleus until NEBD takes place. 

Most Greatwall kinase resides in the nucleus during interphase 
when CdkLCycB activity is low (Burgess et al., 2010; Voets and 
Wolthuis, 2010) suggesting that unphosphorylated Greatwall 
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mainly accumulates in the nucleus. Since active Cdkl:CycB 
enters the nucleus before NEBD, nuclear phosphorylation of Gwl 
becomes enhanced by co-localization of substrate (Gwl) and 
enzyme (Cdkl:CycB). Assuming that Gwl phosphorylation is 
localized to the nucleus before NEBD, but dephosphorylation is 
not, the change in phosphorylated Gwl concentration ([Gwlp]) 
over time is given by: 



d[Gwlp] 
dt 



-k ga -e-[Cdkl : CycB]-[Gwl\—kgj-[B55]-[Gwlp] 



where k ga and k gi are the rate constants for Gwl activation 
(phosphorylation) and inactivation (dephosphorylation) 
respectively. The 's' factor describes the fold increase in rate 
of Gwl phosphorylation as a consequence of colocalization of 
Cdkl:CycB and unphosphorylated Gwl in the nucleus. The value 
of 'e' is proportional to the inverse of the nucleo-cytoplasmic 
ratio. The nucleo-cytoplasmic ratio is around 10% in eukaryotic 
cells (Huber and Gerace, 2007; Joerger and Fersht, 2007), which 
gives a value for 'e' around 10 in interphase. In principle, at 
NEBD the value of 's' should be reduced to one if nuclear 
CdkECycB and Gwl get quickly distributed throughout the 
whole cell. However, neither of them occupies the entire cell 
volume during mitosis, but rather they preferentially accumulate 
around chromosomes and the spindle (Burgess et al., 2010; Voets 
and Wolthuis, 2010). To reflect this inhomogeneous sub-cellular 
localization, we reduce 8 to 25% upon NEBD. To simplify the 
presentation we scale 's' to one in interphase and reduce it to 
0.25 in mitosis after NEBD. 

ENSA is evenly distributed in the nucleus and in the 
cytoplasm, at least in Drosophila (Drummond-Barbosa and 
Spradling, 2004). Because PP2A-B55 substrates localize both 
in the cytoplasm and in the nucleus, we assume that both ENSA 
and PP2A-B55 are evenly distributed throughout the cell. 

We also supplement the model with a downstream mitotic 
substrate (S) that is phosphorylated and dephosphorylated by 
Cdkl and PP2A-B55, respectively. The phosphorylation state of 
this substrate determines whether the cell is in G2 phase 
(unphosphorylated) or in mitosis (phosphorylated). We use the 
phosphorylation state of this substrate (substrate-P or S-P) as a 
proxy for nucleolin phosphorylation that is measured in the 
experiments of Potapova et al. (Potapova et al., 2011). We also 
assume in the model that NEBD takes place when the 
phosphorylation state of S exceeds 70%. 

All the above assumptions are made to simplify the model 
without compromising its main purpose: to understand the 
dynamical processes of normal mitotic entry and mitotic collapse. 

Results 

Temporal dynamics during normal mitotic entry and mitotic 
collapse 

The dynamics of the molecular regulatory system controlling 
mitotic entry in HeLa cells was recently characterized by Potapova 
et al. after synchronization in S phase with double thymidine block 
(Potapova et al., 2011). Most cells completed DNA replication 
around six hours after thymidine wash-out. Further mitotic 
progression of these G2 cells was monitored by means of 
histone-H3 phosphorylation in the presence of a spindle poison 
(nocodazole). In order to investigate the role of Cdkl inhibitory 
phosphorylation, the Weel/Mytl inhibitor (PD0166285) and the 
Cdc25 inhibitor (NSC663284) or their combination was applied 



six hours after thymidine release, i.e. at the beginning of G2 phase. 
Most G2 cells entered into mitosis within 10 hours after thymidine 
release in the absence of any chemical inhibitor. The Cdc25 
inhibitor blocked mitotic entry during this time window, which is 
indicated by the lack of histone-H3 phosphorylation during the 
next four hours after S phase completion. On the other hand, Weel/ 
Mytl inhibition accelerated histone-H3 phosphorylation, and 
mitotic entry were seen within one hour after S phase 
completion. Adding the two (Weel/Mytl and Cdc25) inhibitors 
together also advanced histone-H3 phosphorylation, but only 
transiently because the mitotic state collapsed. 

We have used our model to simulate the time-course 
experiments in figure 5 of Potapova et al. (Potapova et al., 
201 1). Numerical simulations were started with initial conditions 
characteristic for cells in early G2 phase (Weel/Mytl and PP2A- 
B55 active, Cdkl:CycB, Cdc25 and Gwl inactive, and low level 
of CycB etc.; supplementary material Table SI). The initial CycB 
level (0.5 a.u.) was not enough to activate Cdkl in control cells. 
Cdkl:CycB activation (measured by HI kinase activity in the 
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Fig. 2. Time series simulations of mitotic progression and mitotic collapse. 

The time series simulations start from early G2 state with low CycB levels, low 
Cdkl:CycB activity and high PP2A-B55 activity. The temporal patterns of 
mitotic regulators are plotted in four different situations: normal mitotic 
progression (blue), Weel/Mytl (green), Cdc25 (red) and Weel/Mytl and 
Cdc25 double (orange) inhibition. Phosphorylated mitotic substrate 
(A), phosphorylated Weel, Cdc25 and Gwl (B), fraction of active Cdkl 
(C), active PP2A-B55 (D), fraction of inactive (phosphorylated) Cdkl:CycB 
(E) and active (phosphorylated) Ensa-P (F). The model simulations are plotted 
in this way to facilitate the direct comparison to the following experimental 
figures in Potapova et al.: figure 5A,C (A), figure 5B (B), figure 6A (C) and 
figure 5C (E) (Potapova et al., 201 1). Model simulations in panels D and F have 
not been confirmed experimentally yet, therefore they are predictions. 
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experiments) and mitotic substrate phosphorylation required 
three hours of CycB synthesis during G2 phase (Fig. 2A,C). 
Cdc25 inhibition prevented phosphorylation of mitotic substrate 
and Cdkl activation despite increasing CycB levels. Inhibition of 
Weel/Mytl with or without the Cdc25 inhibitor causes premature 
substrate phosphorylation and mitotic entry. Since Weel/Mytl 
inhibited cells can enter into M phase with the CycB level present 
initially (0.5), we simply simulated their mitotic progression with 
constant levels of CycB, which is consistent with the experi- 
mental data. Mitotic substrate phosphorylation (S-P) stays high 
when Weel/Mytl is inhibited which corresponds to a 
stable mitotic state (Fig. 2A). In contrast, mitotic substrate 
phosphorylation is transient in case of the double inhibition 
(Fig. 2A), consistent with the mitotic collapse in the experiments. 
The model simulations have also recaptured the phosphorylation 
patterns of Weel, Cdc25 and Gwl (compare figure 5B of 
Potapova et al. (Potapova et al., 2011) and simulations in 
Fig. 2B). Remember that Weel/Mytl, Cdc25 and Gwl are more 
than just passive readouts of the Cdkl/PP2A-B55 activity ratio, 
because they are also active regulators of the Cdkl:CycB kinase 
and the PP2A-B55 phosphatase. 

During mitotic collapse, the transient phosphorylation of 
Cdkl:CycB substrates results from transient elevation of the 
Cdkl/PP2A-B55 ratio. This transient increase of the kinase/ 
phosphatase ratio is mainly caused by the change in PP2A-B55 
rather than Cdkl activity (Fig. 2C,D). The simulated Cdkl 
activity (Fig. 2C) confirms the experiments in figure 6A of 
Potapova et al. (Potapova et al., 2011), while the change of 
phosphatase activity (Fig. 2D) is a prediction of our model that 
waits for experimental verification. 

To provide a comprehensive picture, the current model is used 
to simulate the temporal pattern of inhibitory phosphorylated 
Cdkl:CycB (pre-MPF, Fig. 2E) and that of phosphorylated 
ENSA (ENSA-P, Fig. 2F). The simulation of pre-MPF agrees 
well with the experimental finding in figure 5C of Potapova et al. 
(Potapova et al., 201 1), while the predicted Ensa dynamics awaits 
further experimental testing. 

In summary, our kinetic model reproduces the temporal 
dynamics that have been carefully observed during normal 
mitotic entry and mitotic collapse. Besides reproducing the 
temporal patterns of experimentally measured mitotic regulators, 
the model also predicts dynamic changes of components that 
have not been followed. However, the numerical simulations 
alone are insufficient to provide a satisfactory mechanism of 
mitotic collapse. To complement the time series simulations and 
to provide a deeper mechanistic insight, we analyse the model 
with the mathematical tool of bifurcation diagrams. 

Explanation of mitotic collapse by dynamical systems theory 
It is informative to analyse the model for mitotic entry using one- 
parameter bifurcation diagrams. A one-parameter bifurcation 
diagram plots stable and unstable attractors (e.g. steady states) of 
a dynamical system as functions of an arbitrary 'parameter'. We 
choose CycB as a bifurcation 'parameter'; although it is not 
constant during the cell cycle, it is the slowest dynamic variable 
of the mitotic control network. The influence of CycB 
(bifurcation parameter) on the dynamical network can be 
illustrated by any of the molecular regulators. We choose 
arbitrarily the generic mitotic substrate (S), which corresponds 
to phosphorylated nucleolin and histone-H3 in the experiments, 
to characterize the state of the system (Fig. 3A). At low CycB 
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Fig. 3. Bifurcation analysis of mitotic progression and mitotic collapse. 

One- and two-parameter bifurcation diagrams are shown for normal mitotic 
progression (A,B), Weel/Mytl inhibition (C,D) and Weel/Mytl and Cdc25 
double inhibition (E,F). One parameter bifurcation diagrams (A,C,E). CycB 
level is used as the bifurcation parameter while phosphorylated mitotic 
substrate indicates the state of the mitotic controls system. The black and the 
red curves represent the bifurcation diagrams before (e=1) and after (e=0.25) 
nuclear envelope breakdown, respectively. Solid lines represent stable while 
dashed lines unstable steady states. Normal mitotic progression (A), Weel/ 
Mytl inhibition (C) and Weel/Mytl and Cdc25 double inhibition (E). The 
cyclin thresholds (labelled by Greek 9) for mitotic substrate phosphorylation 
and dephosphorylation are indicated with subscripts one and two, respectively. 
Between the two corresponding cyclin thresholds the mitotic susbtrate has two 
alternative stable steady states (G2 and M). Two parameter bifurcation 
diagrams (B,D,F). CycB and epsilon (e) are chosen as the two parameters. The 
left black curve records the values of the cyclin thresholds (9 2 values) for the 
mitotic substrate (S) dephosphorylation at different e values. Its value is 0 2 
when s=l and its value is 0 2 ' when s=0.25. The right black curve records the 
values of the cyclin thresholds for the mitotic susbtrate phosphorylation, which 
is 0! when s=l and 0/ when e=0.25. The bistability region is between these 
two curves. The red lines record the time dependent trajectory of the system 
during mitotic progression. During normal mitotic progression (B) the 
trajectory starts from the G2 state in the bistable regime (red dot) and moves to 
the right at constant epsilon (e) value until it crosses the cyclin threshold for 
phosphorylation of the mitotic substrate. As cells enter mitosis, NEBD triggers 
epsilon to decrease so the trajectory drops vertically and ends up the M phase 
state of the bistable regime. After Weel/Mytl inactivation (D), mitotic entry is 
triggered without the requirement of CycB production and the only stable 
steady state is in M phase (red dot). NEBD causes s to decrease which drives 
the system vertically down, but it does not leave the M phase regime. When the 
two inhibitors for Weel/Mytl and Cdc25 are added after S phase (F), the 
system trajectory also drops vertically. The decrease of epsilon (e) causes the 
system to enter the region with unphosphorylated mitotic substrate 
(mitotic collapse). 
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levels, the mitotic substrate is unphosphorylated therefore this 
steady state refers to the G2 phase of the cell cycle. At high CycB 
levels, Cdkl is active and mitotic substrate is fully 
phosphorylated as in M phase (Fig. 3A, black curve). The 
positive circuits of the regulatory network define a bistable 
switch with two different CycB thresholds for mitotic entry (0i) 
and mitotic exit (0 2 ). Between these two cyclin thresholds the 
control system can occupy either of the two stable states (G2 or 
M), and it will depend on the history of the system as to which 
one it settles on. For that reason, the bistable switch is also called 
as a hysteresis switch. The CycB threshold for mitotic entry is 
sensitively dependent on the level of Greatwall-kinase; therefore, 
our model is consistent with the experimental data showing that 
Gwl reduces the Cdkl activity required for NEBD (Hara et al., 
2012). 

The black curve on Fig. 3A characterizes the system with 
unphosphorylated Gwl concentrated in the nucleus. Once the cell 
reaches the 6! cyclin threshold for the G2/M transition, it starts its 
vertical movement to the upper steady state corresponding to M 
phase. During mitotic entry the nuclear envelope breaks down, 
which slows down Gwl phosphorylation because of dispersion of 
Cdkl: CycB and Gwl and weakens PP2A-B55 inhibition within 
the cell. The drop in PP2A-B55 inhibition makes activation of 
Cdkl more difficult, which is manifested in the increase of both 
cyclin thresholds (0^ and 0 2 ' of the red curve). In the model, we 
use the decrease of a single parameter (s) to represent the effects 
of NEBD. Therefore the cellular response to NEBD is better 
illustrated on a two-parameter bifurcation diagram where we plot 
the two cyclin thresholds (the left and the right edges of the 
S-shaped curve) over a range of s values (Fig. 3B). 

On this diagram, the left curve connects the points of mitotic 
exit thresholds (8 2 ) at different values of s. Therefore, the area 
left and below this curve represents G2 phase of the cycle where 
mitotic substrate (S) is unphosphorylated. The right curve 
corresponds to the cyclin thresholds for mitotic entry (0[) and 
the area above and to the right of this curve represents M phase 
where the mitotic substrate is phosphorylated (S-P). Between the 
two curves the system is bistable with low (G2) and high (M) 
mitotic substrate phosphorylation states coexisting. This should 
be imagined as a pleated sheet from a top view where the folds of 
the sheet correspond to the curves on the diagram. 

Early G2 phase cells (CycB=0.5 and £=1) start their mitotic 
journey from the low steady state in the bistable region because 
they have sub-threshold CycB levels. They move to the right by 
CycB synthesis (red horizontal line) and they enter mitosis once 
they pass the 61 curve. Since mitotic entry is accompanied by 
NEBD, which weakens nuclear phosphorylation of Gwl, e gets 
decreased (vertical red line). The system falls back into the 
bistable region, but the mitotic substrate remains phosphorylated 
(S-P) because the 6 2 threshold is not crossed and the mitotic state 
persists in the bistable regime. 

When Weel/Mytl kinases are inhibited both CycB thresholds 
(0! and 6 2 ) are decreased and the bistable region also becomes 
greatly reduced (Fig. 3C). Since the cyclin threshold for mitotic 
entry (61) is lower, cells can now enter into mitosis at a CycB 
level (CycB=0.5) which was not permissive in the absence of the 
inhibitor. Once the phosphorylation of mitotic substrate reaches a 
threshold (0.7 in the model), NEBD takes place. Now both cyclin 
thresholds are increased due to the drop in e, however this 
increase is relatively slight (Fig. 3C, red curve). Since the cyclin 
threshold for mitotic entry (0i) stays below the actual CycB level, 



cells keep on progressing towards highly phosphorylated mitotic 
substrate state and settle there. On the two parameter bifurcation 
diagram the trajectory becomes a vertical line (Fig. 3D), with 
both initial and final states in the highly phosphorylated substrate 
state (S-P), which is M phase of the cell cycle. In summary, 
inactivation of inhibitory kinases advances mitotic entry at low 
CycB level, and mitotic entry is still irreversible. 

When Weel/Mytl and Cdc25 are both inhibited, cyclin 
thresholds are decreased as well. Since the hysteresis effect 
relies on the positive circuits controlling Cdkl activation, the 
bistable region also becomes greatly reduced (Fig. 3E). Similar 
to the Weel/Mytl inhibition case, cells can now enter into 
mitosis at a lower level of CycB (0.5). At NEBD both cyclin 
thresholds are increased (Fig. 3E, red curve), but in this case the 
cyclin threshold for mitotic exit (0 2 ') becomes larger than the 
threshold for mitotic entry (0]). The gap created between the two 
thresholds (0| and 0 2 ') destabilizes the mitotic state and reverts 
the cell back to G2 phase with unphosphorylated mitotic 
subtrates (S), which corresponds to mitotic collapse. 

Again, the cellular fate is better revealed by plotting the time 
dependent trajectory of the system on the 'CycB-s' two- 
parameter bifurcation diagram (Fig. 3F, red line). The same 
initial conditions as control cells (CycB=0.5 and 8=1) define a 
point in the M phase regime. Activation of Cdkl and 
phosphorylation of mitotic substrates accompany M phase entry 
and cause NEBD. NEBD reduces the value of e which pushes the 
system down through the small region of bistability into the G2 
state with mitotic substrates (S) dephosphorylated. The journey 
from the phosphorylated mitotic substrates region into the 
unphosphorylated one corresponds to phenomena of mitotic 
collapse. 

Suppression of mitotic collapse 

The bifurcation diagram on Fig. 3E provides a clear indication 
that mitotic collapse only happens in a restricted window of 
CycB levels. The CycB level must be larger than the threshold 
for mitotic entry (0i is at CycB=0.35 on the black the curve of 
Fig. 3E) in order to bring about mitosis including NEBD in the 
first place. This requirement determines the lower boundary of 
the restricted CycB range. On the other hand, the CycB level 
must be smaller than the threshold for mitotic exit (0 2 ' is at 
CycB=0.55 at on the red curve of Fig. 3E) after NEBD. The 
lower and the upper limits of CycB levels suggest that mitotic 
collapse can take place after Weel/Mytl and Cdc25 double 
inhibition only in a short window of the cell cycle in early G2 
phase. This conclusion is fully confirmed by experimental data in 
figure 6B of Potapova et al. (Potapova et al., 2011). When cells 
were first arrested in mitosis by nocodazole before treatment with 
Weel/Mytl and Cdc25 inhibitors, the mitotic state was stable for 
at least 4 hours. The time series simulations can nicely recapture 
the results of this experiment as well as predict the temporal 
dynamics of unseen mitotic regulators (Fig. 4). Surprisingly, the 
mitotic state is resistant against Cdc25 inhibition only as well 
(Fig. 4; figure 6B of Potapova et al. (Potapova et al., 2011)), 
because in M phase Weel/Mytl kinase are kept inactive by 
Cdkl:CycB. In summary, these results suggest that high CycB 
level suppresses mitotic collapse because the positive circuits of 
Cdkl:CycB activation are strong even if their strength is 
compromised by inhibitors. 

Potapova et al. were also able to suppress mitotic collapse by 
inhibiting PP2A with okadaic acid (OA) (Potapova et al., 2011). 
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Fig. 4. Simulating the responses of mitotic cells to chemical inhibitors. 

(A-F)The model is first allowed to reach a mitotic state, with high CycB level 
(1.3) and low PP2A-B55 activity, and time series simulations are started from 
this steady state. As in Fig. 2, different panels record the dynamics of different 
molecules and four colours are used to represent four different experimental 
conditions. The dynamics of phosphorylated mitotic substrate (A) is directly 
comparable to figure 6B of Potapova et al. (Potapova et al., 201 1). All the other 
panels are model predictions. 

Synchronized HeLa cells at the S/G2 boundary were treated with 
inhibitors of Weel/Mytl and Cdc25 followed by okadaic acid 
(OA) one hour later. Addition of the PP2A inhibitor (OA) was able 
to stabilize the mitotic state and thereby suppress mitotic collapse 
at least for the next 4 hours (figure 6C of Potapova et al. (Potapova 
et al., 2011)). This experiment strongly suggests that mitotic 
collapse is the consequence of misregulated 'activation' of PP2A 
during M phase. 

The time series simulations with the model recapitulate the 
experimental results, as PP2A inhibition stabilises phosphorylation 
of mitotic substrates (supplementary material Fig. SI). In addition, 
the model predicts sustained phosphorylation of Wee 1, Cdc25 and 
Gwl in the absence of PP2A-B55 activity (supplementary material 
Fig. SI). Despite the almost complete phosphorylation of Wee 1 and 
Cdc25, the HI kinase activity is predicted to be low due to the 
inhibition of Cdc25 by NSC663284 (supplementary material Fig. 
SIC). Partial Cdkl dephosphorylation observed in figure 6C of 
Potapova et al. (Potapova et al., 2011) is also captured by our time 
series simulation (supplementary material Fig. SID). 

Discussion 

During progression through the cell cycle, cells undergo multiple 
transitions and qualitative changes to their physiological state. In 



general, there are five transitions during the eukaryotic cell cycle: 
Start (Restriction Point), Gl/S, G2/M, meta/anaphase transitions 
and mitotic exit. It is important that all of these transitions are 
unidirectional and the cell does not revert back to the previous 
phase because that would jeopardise the strict alternation of 
chromosome replication and then segregation. In our view, all the 
eukaryotic cell cycle transitions are made unidirectional 
(irreversible) by underlying bistable switches (Novak et al., 
2007; Verdugo et al., 2013). This proposal has been validated 
experimentally for all the cell cycle transitions of the budding 
yeast cell cycle (Charvin et al., 2010; Cross et al., 2002; Holt 
et al., 2008; Lopez-Aviles et al., 2009) and for the G2/M 
transition in Xenopus oocytes (Pomerening et al., 2003; Sha et al., 
2003). It is an important question whether the design principle of 
bistable cell cycle transitions is an evolutionally conserved 
property of the eukaryotic cell cycle control mechanisms. For this 
reason, we have studied the dynamics of the G2/M transition of 
mammalian cells which was carefully analysed experimentally 
by Potapova et al. (Potapova et al, 201 1). Bistability of the G2/M 
transition could arise from the positive circuits regulating 
CdkliCycB activation through Weel/Mytl and Cdc25. 
Inhibition of both Weel/Mytl and Cdc25 early in G2 phase 
advanced cells into mitosis but progression though M phase was 
aborted. This has been named mitotic collapse (Potapova et al., 
2011). Since inhibition of the enzymes responsible for inhibitory 
phosphorylations weakens bistability, this observation already 
supports the existence of a hysteresis switch. However it remains 
to be answered why cells exit prematurely from mitosis and what 
is the role of phosphatase regulation in this process. 

According to our model, mitotic entry of mammalian cells is 
controlled by a bistable switch controlling Cdkl:CycB activation 
and inactivation of its counter-acting phosphatase, PP2A-B55. In 
G2 phase of the cell cycle, Cdkl is kept inactive by inhibitory 
phosphorylation while PP2A-B55 is active. The accumulation of 
CycB pulls the switch when the threshold for mitotic entry is 
reached at the G2/M transition. Therefore the rise of CycB is a 
trigger pulling signal that turns on the switch by activating Cdkl 
and inhibiting PP2A-B55. During the process of mitotic entry, 
nuclear envelope breakdown (NEBD) creates a cellular 'stress' 
that can jeopardise the G2/M transition. We propose that NEBD 
slows down the phosphorylation of Greatwall-kinase, which is 
critical for effective inhibition of PP2A-B55 and therefore of 
mitotic entry. During mitotic entry, this weakening of PP2A-B55 
inhibition is normally compensated by the positive circuits 
driving Cdkl:CycB activation. If the positive circuits are 
compromised by Weel/Mytl and Cdc25 inhibitors, then PP2A- 
B55 overcomes Cdkl activity and the mitotic state gets 
destabilized prematurely. We are not claiming that the slow- 
down of Gwl phosphorylation by NEBD is the only stress during 
mitotic progression. Another obvious 'stress' during mitotic 
progression is the degradation of Cyclin A by APC/C, since 
Cdkl:Cyclin A is helping Cdkl:CycB to be activated before 
mitosis. Cyclin A degradation is not controlled by the mitotic 
checkpoint and it takes place during prometaphase which implies 
that Cdkl:CycB alone must then maintain the high Cdkl state for 
further mitotic progression. Nevertheless, we cannot conclude 
that mitotic collapse is caused by Cyclin A degradation, as its 
level does not appear to drop during mitotic collapse (figure 5C 
of Potapova et al. (Potapova et al., 2011)). 

Interestingly enough, fission yeast cells with simultaneous 
inactivation of Wee 1 and Cdc25 also undergo occasional mitotic 



Mitotic latching switch 930 



collapses (Sveiczer et al., 1999). Cells of weel s cdc25A double 
mutant show quantized cell cycles with random transitions 
among one, two and three units of cycle time (Sveiczer et al., 
1996). The extended cycle times are the consequences of one or 
two unsuccessful attempts to execute mitosis after which cells 
reset back to early G2 phase and try to enter into M phase again 
later (Sveiczer et al., 1999). The probability of successful mitosis 
increases with cell size which presumably correlates with higher 
level of the Cdcl3 mitotic cyclin. 

We speculate that the requirement of bistability to suppress 
cellular stresses during cell cycle transition is not specific for the 
G2/M transition, but rather a generic phenomenon of the 
eukaryotic cell cycle control. Another example is the meta-to- 
anaphase transition, when cohesins holding sister-chromatids 
together are cleaved. Cohesins are important to provide tension 
against bipolar spindle forces during prometaphase and the lack 
of tension activates error correction and the mitotic checkpoint at 
this stage of mitosis. Yet cohesin cleavage does not reactivate 
error correction and the mitotic checkpoint during anaphase. This 
strongly suggests an underlying bistable switch (He et al., 201 1). 
We propose that cellular 'stresses 1 accompanying eukaryotic cell 
cycle transitions may have been the selective forces necessary for 
underlying bistable switches to have evolved. 

Materials and Methods 

Model simulation and purpose driven parameter estimation 

During the model construction, the molecular network in Fig. 1 was converted into 
a set of ordinary differential equations. To avoid unjustified nonlinearity, the 
model uses mass action kinetics only. The purpose of the current model is to 
understand the quantitative dynamics and the qualitative picture of mitotic entry 
and mitotic collapse. To achieve this goal, we adapt a trial and error method to 
estimate the model parameters. Starting from an initial set of guessed values, we 
adjust the parameters by hand until the model simulations agree with the 
experimental observations (Potapova et al., 2011). Once they agree, we stop 
parameter twiddling. In this regard, our current parameter set is a "working set" 
but by no means optimal, we supply the model code for people who are interested 
in further optimizing the parameters. 

Time series simulations and bifurcation analysis of the model are carried out 
with XPPAUT, a freely available software (http://www.math.pitt.edu/~bard/xpp/ 
xpp.html). 

For the model simulations, we assume that nuclear envelope breaks when 
Cdkl:CycB has phosphorylated the majority (>70%) of the mitotic substrates. We 
reduce 8 to 25% at the time of NEBD. 

We assume that the triggering signal of mitotic entry is the accumulation of 
CycB. Figure 5B of the experimental paper (Potapova et al., 2011) indicates that 
CycB is mainly produced during G2 phase. Furthermore, CycB production seems 
to be blocked after the cell enters mitosis or when the Weel/Mytl inhibitor is 
added. We set the model simulations accordingly. CycB level increases from 0.5 to 
one in absence of Weel/Mytl inhibitor, and CycB level is kept at 0.5 in the 
presence of Weel/Mytl inhibitor. 

The experimental paper has reported measurements of HI kinase activity and of 
pre-MPF. In our model, we represent HI kinase activity by the ratio between the 
active, unphosphorylated Cdkl:CycB and the total CycB. Pre-MPF is represented 
by the ratio between phosphorylated Cdkl:CycB and the total CycB. 
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